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Abstract The boundary layer analysis of a 2D forced convection flow along a persistent moving 
horizontal needle in electrically conducting magnetohydrodynamic dissipative ferrofluid is investi¬ 
gated. The energy equation is constructed with the Joule heating, variable heat source/sink and dis¬ 
sipation effects. To check the variation in the boundary layer behaviour, we considered the two 
ferrofluids namely, Fe 3 0 4 -methanol and Fe 3 0 4 -water. The reduced system of governing PDEs 
are solved by employing the R-K process. Computational results of the flow and energy transport 
are interpreted with the assistance of tabular and graphical illustrations. Increasing the needle size 
significantly reduces the flow and thermal fields of both nanofluids. In particular, thermal and 
velocity fields of Fe 3 0 4 -methanol nanofluid is highly depreciated when equated with the Fe 3 0 4 - 
water nanofluid. 

© 2017 Faculty of Engineering, Alexandria University. Production and hosting by Elsevier B.V. This is an 
open access article under the CC BY-NC-ND license (http://creativecommons.Org/licenses/by-nc-nd/4.0/). 


1. Introduction 

Recent days, the momentum and heat transfer over an inserted 
solid object has attracted the minds of several investigators 
because of their significant applications in several industries 
and engineering areas. Few of them are a hot wire anemometer 
or covered thermocouple are used to quantify the airflow 
velocity, used as heat remover from micro scale equipment. 
In view of these applications, Chamkha and Issa [1] elaborated 
the impact of heat source/sink and thermophoresis on 2D 
magnetohydrodynamic flow past a semi-infinite porous flat 


Corresponding author. 

E-mail address: maths.sulochana@gug.ac.in (C. Sulochana). 

Peer review under responsibility of Faculty of Engineering, Alexandria 
University. 


surface with the non-uniform magnetic field. A similar kind 
of investigation was carried by Silva and Gosselin [2], They 
incorporated an interesting flow geometry for the higher heat 
transport rates. Further, non-uniform flow of micropolar fluid 
through a circular cylinder was numerically analysed using 
Keller-box method by Ali et al. [3]. Later on, several authors 
(Refs. [4-8]) proposed the mathematical models to understand 
the momentum and energy transfer over various body shapes 
such as a parallel plate, vertical cone, circular cylinder, sphere, 
wedge, and channels. 

Geometrically, the non-uniform thickness is one of the 
important characteristic features of body objects. An occur¬ 
rence of these studies can be found in space science and indus¬ 
tries. Historically, the investigation of flow past an object with 
variable thickness was initiated by Lee [9] and he discussed the 
momentum and thermal transport behavior over a needle. 


http://dx.doi.Org/10.1016/j.aej.2017.08.020 

1110-0168 © 2017 Faculty of Engineering, Alexandria University. Production and hosting by Elsevier B.V. 
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Nomenclature 



U , V 

velocity components in x and y — directions 

Greek letters 

x, r 

cylindrical coordinates, in axial and radial direc- 


surface shear stress 


tions 

e 

dimensionless temperature 

U 0 

composite velocity 

<t> 

solid volume fraction 

Mo, Mqq 

velocity of the needle and velocity of the main 

Pnf 

density of the nanofluid 


stream 

ftnf 

dynamic viscosity of the nanofluid 

To, ^ 

temperature at the surface and ambient 

Gnf 

electrical conductivity of the nanofluid 

c 

thickness of the needle 

Pf 

density of the base fluid 

R(x) 

surface shape of the thin needle 

Ps 

density of the nanoparticle 

k f 

thermal conductivity of the base fluid 

V 

kinematic viscosity 

ks 

thermal conductivity of the solid particles 

( PCp)nf 

specific heat capacity of the nanofluid at constant 

Pr 

Prandtl number 

pressure 

M 

magnetic field parameter 

(pCp)f 

heat capacity of the base fluid at constant pressure 

Ec 

Eckert number 

(. pC P ), 

heat capacity of the nanoparticle 

A * 

space dependent heat source/sink 

+ 

stream function 

B* 

temperature dependent heat source/sink 

n 

similarity variable 

Cf 

Nu x 

skin friction coefficient, 
local Nusselt number 

y 

velocity ratio parameter 


Then after, few authors (Narain and Oberoi [10]) extended the 
Lee [9] work by considering forced convection fluid on a thin 
needle. Later, Lee et al. [11] examined the mixed convection 
flow over two different geometries with uniform heat flux. Fur¬ 
ther, the related type of study was reported by Ahmad et al. 
[12] for two cases such as assisting flow and opposing flow. 
Since then, many investigators (Sulochana and Sandeep [13], 
Khan et al. [14,15] Babu and Sandeep [16]) worked on fluid 
flow over a thin needle (i.e. non-uniform thickness), slendering 
sheet, paraboloid of revolution, stretching sheet. 

Generally, conventional fluids used in factories such as 
water, kerosene, transformer oil cannot meet prerequisites of 
effective thermal conductivity as needed by the industries. 
The use of solid particles of the range 1-100 nm in conven¬ 
tional fluids improves the thermal conductivity drastically. 
These minute solid molecules are called nanoparticles. And 
use of nanoparticles in a conventional fluid is called nanofluid. 
Scientifically, use of nanomaterial in conventional fluid 
upgrades the thermal conductivity to a larger extent, at the 
very first this concept was introduced by Choi [17]. It is well 
known that two basic models to describe nanofluid character¬ 
istics are Buongiorno [18] and Tiwari and Das [19], To describe 
the nanofluid behavior Buongiorno [18] introduced a very 
strange model with incorporating thermophoresis and Brown¬ 
ian motions. The author found that thermophoresis and Brow¬ 
nian motions are more effective among all used slip 
mechanisms. Some of the authors (Zulfiqar et al. [20], Hassani 
et al. [21]) used this model to describe the nanofluid behavior. 
One more model developed by Tiwari and Das [19] is also 
helpful in describing nanofluid behavior. Again, this model is 
also being used by Soid et al. [22]. Very recently, Sheik- 
holeslami and Zeeshan [23] numerically examined the impact 
magnetic field and constant heat flux on the flow of water 
based CuO nanofluid by using CVFEM. They incorporated 


KKL model and found that improving buoyancy forces corre¬ 
spondingly decay heat transfer rate. 

Moreover, the impact of radiation on MHD flow of vis¬ 
cous, incompressible fluid through a permeable stretching/ 
shrinking surface with nanoparticles was interrogated by 
Ibrahim and Shankar [24]. Sandeep et al. [25] presented a 
comparative analysis of three different kinds of fluids (i.e. Jef¬ 
frey, Maxwell and Oldroyd-B) over a stretching surface with 
considering the suction/injection effects. They found that 
Oldroyd-B nanofluid is largely reactive in improving Nusselt 
and Sherwood numbers as compared with Jeffrey and Max¬ 
well nanofluids. Wahed et al. [26] elaborated 2D MHD flow 
of incompressible nanofluid past a continuous moving sur¬ 
face. For this study, they have taken various effects into 
accounts such as variable thickness. Brownian motion and 
heat generation. A very good number of articles are available 
on the topic nanofluid flow over various geometries with dif¬ 
ferent nanoparticles, few of them are [27-34]. Scientifically, 
the fluids with magnetic nanoparticles are called ferrofluids. 
These fluids exhibit more thermal conductivity as compared 
with other ordinary nanofluids. Sandeep [35] examined the 
impact of an aligned magnetic field on the flow magnetic- 
nanofluids (ferrofluids) over a stretched surface. Recently, 
Sheikholeslami [36] discussed the impact of Coulomb forces 
on the flow of water based Ferro fluids for the thermal 
enhancement. 

In this analysis, we incorporated the Tiwari and Das model 
[19] to investigate the boundary layer analysis of a 2D forced 
convection flow along a persistent moving horizontal needle 
in electrically conducting magnetohydrodynamic dissipative 
ferrofluid. To check the variation in the boundary layer behav¬ 
ior, we considered the two ferrofluids namely, Fe 3 0 4 -methanol 
and Fe 3 0 4 -water. The reduced system of governing PDEs are 
solved by employing the R-K based shooting process. 
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Fig. 2 Velocity behavior with Magnetic field parameter M. 

2. Mathematical formulation 

We ponder an electrically conducting, viscous, incompressible 
2D dissipative ferrofluid flow over horizontally moving needle 
as displayed in Fig. 1. Here (x, r) represent the axial and radial 
coordinates in cylindrical form, respectively, where c is the 
thickness of the needle. 

The size of the needle is thin, as its thickness is not more 
than the boundary layer formed on it. Also, the velocity of 
the horizontally moving thin needle is u„. The flow governing 
equations using Tiwari and Das [19] nanofluid model in cylin- 


drical coordinates are: 


M* 

+ H, = 

= 0 , 

(1) 

PnA U 

(fix) + v( 

fir)} = Pnf(Mr){r(fir)} r ~ OnfB\u, 

(2) 

c pc P : 

)nfM T x) 

l + v(r,)} = yi//){r(T)},. + f 




+ <T„fB^U 2 + f-tffir) , 

(3) 



0 0.5 1 1.5 2 2.5 3 

V 

Fig. 3 Temperature behavior with magnetic field parameter M. 

the appropriate boundary restrictions are: 
u = i/ 0 , v = 0, T = T 0 at r = R(x) "1 

_ f ? v v 

u -> Uoo, T x as ;• -> oo j 

where u and v are the velocity components along axial and 
radial coordinates. And R(x) represents the surface shape of 
the thin needle. The nanofluid parameters are given as 

Krf (k s +2kf)-2<f>(kf-k s ) a n f f ] \ 3(<r-l)(ft \ 'j 

kf {k s +2kf)+<j)(kf-k s ) ’ Gf y (<T+2)-(tT-l)</>y ’ 

a =%’Pnf- PA( l - 4>)+<t>{P s /Pf)},Pnf= Pfi X - 4>)~ 2 ' 5 , f ’ 

( P C p)nf= (P C p)f{A -<t>) + <KP C p)J(P C p)f} } 

(5) 

to reduce the Eqs. (1) (3) to the dimensionless form, the sim¬ 
ilarity transformations used are: 
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7 

Fig. 4 Velocity behavior with volume fraction parameter cf>. 
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Fig. 6 Temperature behavior with A*. 



Fig. 5 Temperature behavior with volume fraction parameter 4>. 


V 

Fig. 7 Temperature behavior with B". 


<A = Vfxf{ri),6= (T - r oo )/(r 0 - T oo),ri = U 0 r 2 /v f x, 'j 
u = (l/r)tl/ r ,v= (-l/r)\l/ x ,R(x) = (v f xc/U)' /2 , > (6) 

q'" = (, kfU 0 (T 0 - Too)/x\f){A*f + B*9}, j 

By setting t; = c in Eq. (6) describes the size and shape of 
the thin needle r = R(x). 

Eq. (1) is automatically satisfied, further the Eqs. (2) and (3) 
reduced to: 


A - (1 — 4>) 5/2 ,B = 1 -(/> + (t>{pJp f ),K = k nf /k f , 
E=\-<t>+4>{(pC p )J{ P C p ) f } 1 (9) 

transformed boundary conditions are: 

A?) = ( 7 / 2 )c,/(c) = y/2, 9(c) =1 "1 

/('/) (1 - 7/2), 9{n) -» y/2 as >1 -> 00 J ’ 


(2Af"'ri + 2zi/") + B(ff" - Mf) = 0 (7) 

2 K{B" n + 9') + EPr0 + (. A*f + B?9) + PrE c (2 n (f"f 

+ M(f') 2 ) = 0, (8) 


where “prime” denotes derivative w.r.t. t; and the dimension¬ 
less parameters are expressed as, 

M = a f Bl/2U 0 p f , Pr = p f C p /k f , E c 

= (2t/o)~ /Cp(To — 74), y = uo/Uo, Uo = Uo + Uoo (11) 


where 
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Fig. 8 Temperature behavior with viscous dissipation parameter 
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Fig. 10 Temperature behavior with needle thickness parameter 
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Fig. 9 Velocity behavior with needle thickness parameter c. 

The physical quantities of practical interest are Cf and Nu x 
which are given by 

C f = T w /p f U 2 0 , Nu x = xq„/kf{T 0 - T x ) (12) 

where x„ and q w , are given as: 
tV = l “nA‘r)r=c = Ul/VfX)f{c) , 

q„, = - k nJ {Tf )- r=f = -k„f{2cU 0 {T 0 - r oo )/v / x}0'(c), 

Using Eq. (6), the above equation reduces as 
Re\ /2 C f = 4c 1/2 4/"(c), Re~ x l/2 Nu x = -2 c l/2 K6'(c) (14) 

where Re x = U 0 x/Vf. 



3. Discussion of the results 

Here we discuss the flow and energy transport behavior of two 
distinct ferrofluids under the influence of sundry physical 
parameters. The transformed dimensionless flow governing 
Eqs. (7) and (8) with associated corresponding boundary 
restrictions (10) are numerically solved by employing R-K 
and Newton methods. The impact of pertinent parameters 
viz. M, cf >, A *, B *, E c and c on flow and energy transport behav¬ 
ior of the fluid, with Nu x and Cf are bestowed with the aid of 
graphs and tabular values. Certain parameter values are kept 
constant i.e. cf> = 0.1; M = 1; rj = 1; A* = 1: B* = 1; 
2=1; c = 0.1; Pr = 6.587 (water) and Pr = 7.38 (Methanol) 
in the complete study, excluding the changes are mentioned 
in Figs. 2-9 and Tables 2 and 3. 

The impact of ascending values of M on temperature and 
velocity fields are portrayed in Figs. 2 and 3. We detect a 
decaying in velocity fields with M. Physically, this may happen 
because of the generation of Lorentz force (i.e. opposing force 
to the flow of ferrofluids). The existence of these force may 
lead to decline the fluid velocity and enlarges the thermal 
boundary layer thickness. Figs. 4 and 5 are depicted to observe 
the action of volume fraction parameter cf> on velocity and tem¬ 
perature fields. From these figures, we noticed that the infla¬ 
tion in volume fraction parameter cf> curtails the velocity 
fields but it inflates the thermal fields. Generally, as cf> improves 
the viscosity of the fluid and then highly viscous fluid moves 
with lesser velocity. Hence we observe a declination in fluid 
velocity and enhancement in temperature fields. It is worth 
to mention that the water- Fe^On nanofluid has higher fluid 
temperature as compared with methanol- Fe^Oa, nanofluids. 

The nature of the temperature field for various values of A* 
and B * are displayed in Figs. 6 and 7. We may observe that the 
increasing values of A* and B* improves the fluid temperature. 
Generally, positive values of A* and B * produces heat and neg¬ 
ative values of A* and B* absorbs heat. Hence it is found an 
improvement in fluid temperature. It is fascinating to write 
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Table 1 Thermo physical properties of water, methanol and 
Fe 3 0 4 . 


Parameters 

Water 

Methanol 

Fe 3 0 4 

p (Kg/m 3 ) 

992 

792 

5180 

C P (J/Kg K) 

4076 

2545 

670 

k (W/mK) 

0.605 

0.2035 

9.7 

a (S/m) 

0.005 

0.5 x 10~ 6 

0.74 x 10 6 

Pr 

6.587 

7.38 

- 


that the thermal boundary thickness of water based ferrofluids 
is slightly higher than methanol based ferrofluids. Fig. 8 illus¬ 
trates the impact of Ec on temperature fields. The figure 


reveals that the rising value of Ec improves the fluid velocity, 
this leads to increase in the fluid velocity near the surface 
and as a result friction heat will be produced. Because of this, 
we have seen enhancement in temperature field. Figs. 9 and 10 
are constructed to study the impact of c on flow analysis. It is 
clearly seen that the improving values of c decelerates the pro¬ 
files of the fluid flow. 

Table 1 is drawn to display the thermophysical properties 
of base fluids and nanoparticles. Table 2 and 3 drawn to shed 
light on the response of Nu x and Cf for various pertinent 
parameters for water- Fe 3 0 4 ferrofluid and methanol- Fe 3 0 4 
ferrofluid, respectively. From above table, we found that, skin 
friction of water based ferrofluids are decelerated more than 
that of methanol based ferrofluids. Heat transfer rates can be 
controlled by A* and B*. The volume fraction of nanoparticles 
can improve the local Nusselt number effectively in Fe 3 0 4 - 


Table 2 

Numerical values of C/ and Nu x for Fe 3 0 4 -water nanofluid for various 

flow parameter values. 


M 

(j) A* B* E c 

c 

Fe 3 0 4 -water 





f (c) 

(c) 

1 



-1.292476 

0.852388 

2 



-1.536943 

0.805100 

3 



-1.738413 

0.766528 


0.1 


-1.292476 

0.852388 


0.2 


-1.827584 

1.043389 


0.3 


-3.897519 

1.596699 


1 


-1.292476 

0.852388 


2 


-1.292476 

0.757849 


3 


-1.292476 

0.663310 


1 


-1.292476 

0.852388 


2 


-1.292476 

0.593761 


3 


-1.292476 

0.191463 


0.1 


-1.292476 

0.852388 


0.5 


-1.292476 

0.439325 


0.9 


-1.292476 

0.026281 



0.1 

-1.292476 

0.852388 



0.3 

-2.297010 

1.819543 



0.5 

-3.041710 

2.794746 



Table 3 Numerical values of Cj■ and Nu x for i*e 3 0 4 -methanol nanofluid for various flow parameter values. 

M (j) A* B* E c c 

Fe 3 0 4 -methanol 



/"(e) 

-O'(e) 

i 

-1.323960 

0.955844 

2 

-1.582287 

0.900177 

3 

-1.794326 

0.855024 

0.1 

-1.323960 

0.955844 

0.2 

-1.900611 

1.218736 

0.3 

-4.118343 

1.983086 

1 

-1.323960 

0.955844 

2 

-1.323960 

0.865306 

3 

-1.323960 

0.774793 

1 

-1.323960 

0.955844 

2 

-1.323960 

0.726384 

3 

-1.323960 

0.399975 

0.1 

-1.323960 

0.955844 

0.5 

-1.323960 

0.493412 

0.9 

-1.323960 

0.030980 

0.1 

-1.323960 

0.955844 

0.3 

-2.356545 

2.064814 

0.5 

-3.125189 

3.198403 
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Table 4 

Validation of the numerical technique for Nu x . 




M 

RKN 

RKF 

RKS 

BVP4C 

BVP5C 

1 

0.955844 

0.9558444521 

0.9558444521 

0.9558444521 

0.9558444521 

2 

0.900177 

0.9001772133 

0.9001772134 

0.9001772134 

0.9001772134 

3 

0.855024 

0.8550247642 

0.8550247641 

0.8550247641 

0.8550247642 

4 

0.803241 

0.8032416430 

0.8032416430 

0.8032416432 

0.8032416432 


methanol ferrofluid as compared with the Fe 3 0 4 -water fer- 
rofluid. Eckert number Ec has a tendency to decline the local 
Nusselt number. Table 4 shows the validation of the computa¬ 
tional results by comparing with the other techniques. 

4. Conclusions 

The 2D boundary layer flow of water and methanol based fer- 
rofluids with non-uniform heat generation/absorption and vis¬ 
cous dissipation is numerically analysed. The electrical 
conductivity of ferrofluid is taken into consideration. A few 
important findings of the present study are listed below. 

• A significant hike in heat transfer performance can be deliv¬ 
ered in methanol based ferrofluid when compared to water 
based ferrofluid. 

• Thermal boundary layer thickness of water-Fe 3 0 4 ferrofluid 
is higher than the methanol-Fe 3 0 4 ferrofluid. 

• Water based Fe 3 0 4 ferrofluid moves faster than the metha¬ 
nol based Fe 3 0 4 nanofluid. 

• The parameters A* and B* acts as a heat transfer rate con¬ 
trolling parameters. 

• The parameter cj) improves the local Nusselt number of 
Fe 3 0 4 -methanol ferrofluid when compared with Fe 3 0 4 - 
water ferrofluid. 
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